Axion-Like Particles as Ultra High Energy Cosmic Rays? 
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If Ultra High Energy Cosmic Rays (UHECRs) with E > 4 x 10 19 eV originate from BL Lacertae 
at cosmological distances as suggested by recent studies, the absence of the GZK cutoff can not be 
reconciled with Standard-Model particle properties. Axions would escape the GZK cutoff, but even 
the coherent conversion and back-conversion between photons and axions in large-scale magnetic 
fields is not enough to produce the required flux. However, one may construct models of other 
novel (pseudo)scalar neutral particles with properties that would allow for sufficient rates of particle 
production in the source and shower production in the atmosphere to explain the observations. As 
an explicit example for such particles we consider SUSY models with light sgoldstinos. 



I. INTRODUCTION 

Ultra High Energy Cosmic Rays (UHECRs) with en- 
ergies above the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off 0] were detected in all relevant experiments @-[|, 
suggesting that these particles can not originate at cos- 
mological distances. On the other hand, there are no 
apparent nearby sources in their arrival direction. There- 
fore, something fundamental appears to be missing in our 
understanding of the sources, nature, or propagation of 
UHECRs. 

The small-scale clustering of UHECR events sug- 
gests that the sources are point-like on cosmological 
scales 0. Several astrophysical sources were suggested 
based on the coincidence of the arrival directions of 
some of the highest-energy events with certain astro- 
physical objects ||. For example, a correlation be- 
tween compact radio quasars and UHECRs was sug- 
gested in ||,|ri],[l3| , although other authors found them 
to be insignificant jlO|,fL2) . Recently, a statistically signif- 
icant correlation, at the level of chance coincidence be- 
low 10 -5 , was found with the most powerful BL Lacertae, 
i.e. quasars with beams pointed in our direction The 
identified sources are at z > 0.1, far exceeding the GZK 
distance of -Rgzk ~ 50 Mpc, so that the primary UHE 
particles can not be protons. The photon attenuation 
length for energies around 10 20 eV is of order the GZK 
cutoff distance, primarily due to the extragalactic radio 
backgrounds. While the limiting magnitude of the ra- 
dio backgrounds necessary to absorb UHE photons can 
be determined only by numerical propagation codes [ fL5] , 
one can even now conclude that UHECRs with energies 
around 10 20 eV are very unlikely to be photons. 

The only Standard-Model particles which can reach 
our Galaxy without significant loss of energy are neu- 
trinos. Two different scenarios involving UHE neutri- 
nos have been proposed. In the first, neutrinos pro- 
duce nucleons and photons via resonant Z-production 
with relic neutrinos clustered within about 50 Mpc from 
the Earth, giving rise to angular correlations with high- 
redshift sources [EgI. However, for the interaction rates to 



be sufficiently high, this scenario requires enormous neu- 
trino fluxes and an extreme clustering of relic neutrinos 
with masses in the eV range |L7fl . The second neutrino 
scenario invokes increased high-energy neutrino-nucleon 
cross sections. This could be caused by the exchange of 
Kaluza-Klein graviton modes in the context of extra di- 
mensions [fl8[ or by an exponential increase of the number 
of degrees of freedom in the context of string theory [jl9| . 

Another possibility to avoid the GZK cutoff is a small 
violation of Lorentz-invariance, a hypothesis which can 
not be tested in terrestrial experiments |po| , pT| . 

The GZK cutoff can be avoided also if the UHECRs 
consist of certain new particles. One possibility is a new 
stable massive hadron with a mass around 2-3 GeV [^2j , 
shifting the GZK bound to higher energies E > 10 21 eV 
into a range where no UHECR event has yet been found. 
However, it now appears that these exotic hadrons are 
excluded by laboratory experiments j2j| . 

Therefore, if the UHECRs indeed originate from point 
sources at cosmological distances one is running danger- 
ously short of plausible explanations for how this radi- 
ation can reach us. This perhaps desperate situation 
motivates us to consider other options for new particles 
which can traverse the universe unimpeded at high ener- 
gies. Specifically, we consider the possibility of axion-like 
particles, i.e. electrically neutral (pseudo)scalar particles 
X with a relatively small mass Mx < 10 MeV. 

Such particles must fulfill several requirements to be 
candidates for UHECRs. They must live long enough 
to reach us from a cosmological distance. They must not 
lose too much energy in interactions with the CMBR and 
other background radiations or in extragalactic magnetic 
fields. They must interact sufficiently strongly in or near 
our Galaxy or in the Earth's atmosphere to produce the 
observed UHE events. Finally, their interactions must 
allow for the production of a significant flux at the source. 

We will first consider proper axions and find that they 
seem to be excluded as UHECRs. We then turn to more 
general particles and study their necessary properties to 
fulfill the above requirements. As an explicit example we 
study light sgoldstinos. 
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II. PROPER AXIONS 

Proper axions arise from the Peccei-Quinn mechanism 
to solve the strong CP problem. As such their properties 
are governed by one main parameter, the Peccei-Quinn 
scale or axion decay constant f a ; astrophysical limits im- 
ply fa ^ 10 10 GeV. Axions mix with neutral pions so that 
their mass and interaction strength are roughly those of 
7T°, reduced by f n /f a with f n w 93 MeV the pion decay 
constant. It is easy to see that axions live long enough 
and interact weakly enough with the CMBR to traverse 
cosmological distances unimpeded. By the same token, 
their interaction strength is far too weak to imagine their 
efficient production at the source or their efficient detec- 
tion in the Earth's atmosphere. 

It is less obvious, however, if they could not be pro- 
duced in sufficient numbers by their coherent conversion 
7 — > a in large-scale magnetic fields in the source region, 
and then re-appear as photons in the galaxy by the in- 
verse process. Put another way, one might imagine the 
UHECRs to be photons which traverse the universe in 
the guise of axions. 

The conversion between axions and photons in a large- 
scale magnetic field is essentially a particle oscillation 
phenomenon 24 1. The diagonal elements of the mixing 
matrix involve to 2 and the square of the "photon effec- 
tive mass" within the given medium, the off-diagonal el- 
ement, which induces the mixing, is 2g ai BE where B is 
the transverse magnetic field and E the particle energy. 

The oscillation length £ osc corresponds to the momen- 
tum difference between axions and photons of the given 
energy, m our case E w I0 20 eV. Noting that the effective 
photon mass is much smaller than m a , the momentum 
difference is governed by the axion mass alone so that 
£ osc = ^E/m 2 a = 8.1 kpc (£yi0 20 eV) (meV/m a ) 2 . On 
the other hand, the coherence length ^ of the galactic 
magnetic field is probably less than I kpc. A significant 
conversion rate requires ^ osc < £b, i-e. m a larger than a 
few meV, not in contradiction with current limits. 

The effective mixing angle between axions and photons 
in a magnetic field is given by \ tan(2#) = g aj BE/m 2 = 
0.2( ffa7 /10~ 10 GeV _1 )(i;/10 20 eV)(meV/m a ) 2 . The as- 
trophysical limit on the axion-photon coupling is g ai < 
1(T 10 GeV" 1 so that for m a not much smaller than 

1 meV the mixing angle becomes large. Put another way, 
for g ai near its limit and m a around 1 meV the transition 
rate in the galaxy is not ridiculously small. 

The numbers are much worse for proper axions where 
g ai and m a are related by g ai w a/{2irf a ) and m a w 
mnfn/fa- In this case the mixing angle becomes 8 sa 

2 x lCT 4 (£;/T0 20 eV) (meV/m a ) « 1. Allowing the ax- 
ion mass to be very small, the mixing angle could be- 
come reasonably large. On the other hand, the oscilla- 
tion length then becomes very much larger than is. The 
transition probability is P(a — > 7) = {g^Bls /2) 2 ~ 

3 x 1(T 8 (B/liG) 2 (^s/kpc) 2 (I0 10 GeV// a ) 2 « 1. A sim- 
ilar estimate applies to the source region where the mag- 



netic field could be stronger, but the correlation length 
would be smaller. Therefore, the combined probability 
to produce axions at the source and to convert them into 
photons in the Galactic magnetic field is tiny, perhaps as 
small as P ~ 1CP 16 . Therefore, if one adjusts the pro- 
ton flux from the sources to the observed flux below the 
GZK cutoff, then no significant axion flux will be pro- 
duced above the cutoff. Therefore, proper axions do not 
work for this scenario. 

Even if one dials g aj and m a independently, the num- 
bers look discouraging as one would need a huge rate 
of UHE photon production in the source to compensate 
for small transition rates both in the source and galaxy 
and one would need parameter values near their exclusion 
limits. 



III. GENERIC AXION-LIKE PARTICLES 

Since proper axions are apparently not able to explain 
the UHECR phenomenon, we next turn to a more ex- 
otic new scalar X; a similar analysis for pseudoscalars is 
straightforward. The new particle is assumed to couple 
to gluons and photons via nonrenormalizable interactions 
of the form 



(1) 



Only these two interactions will be important, so we as- 
sume that the coupling to other Standard-Model particles 
are suppressed because, say, they proceed through loops 
or are proportional to small Yukawa constants. 

If Mx < 2m,r = 270 MeV, the dominant decay mode 
is into two photons: 
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because the direct coupling to electrons is suppressed by 
assumption. If this light particle has the energy Ex it 
propagates through the Universe without decay if 
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where V x is essentially identical with the two-photon de- 
cay rate Eq. (Q). Therefore, we need to require 



< 1.6 x lCT 11 GeV" 



E x Z lOMeV y 
10 20 cV v, M x ) [ ' 



if these particles are supposed to reach us from cosmo- 
logical distances. 



*The axion-photon coupling of the previous section was 
based on the normalization C a -y = (g a -, /4)aFF — g a -,aE ■ B. 
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Propagating through the Universe, the light scalar X 
may also disappear by interactions with the CMBR. For 
E x « 10 20 eV, the CM energy is E cm w {2E x uo) 1/2 « 
350 MeV, where luq w 6 x 10~ 4 eV is the average en- 
ergy of relic photons. Pairs of light charged particles A ± 
are produced with the cross section a(X-f — * A + A~) = 
ag 2 /16. With a relic photon number density of about 
400 cm -3 the requirement Rx-y-^A+A- > ^Universe gives 
g 7 < 1 GeV -1 . Similar estimates apply to other possible 
processes like -V7CMB — ► 7 71 " - Therefore, the tiny pho- 
ton coupling required by Eq. (^) guarantees the absence 
of a GZK cutoff for the X particles. 

Both the production of X particles at the source and 
their interaction in the atmosphere require rather large 
cross sections, comparable to strong ones. For X parti- 
cles with the characteristic energy scale g~ l this is pos- 
sible only if the CM energy in the system is close to this 
scale, but not significantly higher so that the effective in- 
teractions (|l|) are still meaningful. Typical CM energies 
of UHECR interactions with nucleons are E cm ss 100- 
300 TeV. We can estimate the interaction cross section 
with nucleons at such energies as 



ox = os 



ax 



The suppression factor 



ax _ (E cm g g y 
a s 



4ira s 



(5) 



(6) 



should not be very small. 

We next turn to the X mean free path (mfp) ix in the 
Earth's atmosphere. Since our particle exhibits strong 
interactions we estimate ix by analogy with the proton 
mfp £ p as ix = ip {a s /a g ). To initiate an atmospheric 
shower, X should have a relatively small mfp. Assuming 
ix < 10 i p and using Eq. (||) and a s — 0.1 we estimate 



g g > 1.1 x 10~ 6 GeV" 1 



10 20 eV 



(7) 



The inequalities (||) and (fjj) determine the g g range suit- 
able for explaining the UHECRs above the GZK cutoff. 

How are the .V-particles produced at an astrophysical 
source like a quasar? If our estimate for the cross sec- 
tion Eq. (||) is valid, UHE X particles will be efficiently 
produced in the high-energy tail of the proton spectra 
by proton-proton collisions while their production at low 
energies will be negligible. Therefore, we can expect that 
the proton flux from the source at low energies will con- 
tinue with the same slope at high energies due to the X 
component. Only part of the initial proton energy will 
be transfered to the X particles; probably they will be 
produced on the peak of the gluon distribution function 
with E w 0.1i?p. However, once produced they will es- 
cape more easily from the source compared with protons 
precisely because their cross section is smaller. 



Many bounds on axion-like particles arise from cosmol- 
ogy, astrophysics and laboratory measurements p5| , p6| . 
Still, there remain regions in parameter space where 
X particles can explain UHECRs without contradicting 
these limits. In Fig. [I] we present the experimentally 
allowed regions in the space (g 7 ,Mx) where the inequal- 
ity (|4|) is satisfied. In each concrete model one can evalu- 
ate the effective coupling constant g~ f which has to belong 
to the allowed regions shown in Fig. [I]. Since generally 
the interaction with gluons leads at higher order to an ef- 
fective interaction with photons, the inequality (R) may 
shrink the allowed regions in Fig. ^in concrete models. 




FIG. 1. The allowed region for the parameters (Mx, J T ) 
are shaded in grey. The region traced by the long-dashed 
line is ruled out by the helium-burning life-time of horizon- 
tal-branch stars fe4|. The region surrounded by a thin solid 
line is ruled out by SN 1987A. The region confined between 
short-dashed lines is ruled out by the photon background and 
the CMBR EJ. Below the thick solid line the inequality (|) 
is valid. 

From the general case one can see that constraints on 
the X particle interactions favor a strong coupling to glu- 
ons and a tiny one to photons. Hence the first extreme 
example is a light scalar X which interacts at tree level 
only with gluons according to Eq. ([!]); a similar analy- 
sis applies to a light pseudoscalar. The interaction with 
all other SM particles arises at higher order. In particu- 
lar, because the gluonic operator creates mesonic fields, 
the interaction X77 emerges with a coupling constant re- 
specting the hierarchy g^j 1 g g ~ a/(Air) <~ 10~ 3 . In view 
of this relationship the inequality (pn allows only the re- 
gion of parameter space which corresponds to the up- 
per shaded region in Fig. [IJ. Unfortunately, this allowed 
region corresponds to a fairly small g^ 1 ~ 0.1-5 TeV. 
Therefore, our nonrenormalizable model for V-baryon 
scattering in the atmosphere becomes invalid because it 
should proceed at 100 TeV in the CM frame. 

This example shows that the lowest region in Fig. |l| 
is unphysical, because the condition (Q) requires the hi- 
erarchy g 1 jg g ~ 10~ 10 , which is impossible due to loop 
contributions. The Mx ~ MeV region in Fig. [j] can still 
exist in models with a hierarchy between photon and 
gluon couplings, but this requires a two order of mag- 
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nitude fine-tuning for the ratio g~//g g down to values of 
order 10~ 5 . 

The other possibility is that the couplings to photons 
and to gluons are of the same order. In this case only the 
upper region in Fig. [I] is interesting because the gluon 
coupling should not be too small from Eq. (R). We now 
turn to an explicit example for a model which does not 
need any fine tuning of the couplings g 1 and g g . 



IV. LIGHT SGOLDSTINOS 

As an example of a realistic model for X parti- 
cles we consider the supersymmetric extension of the 
SM with a light scalar and/or pseudoscalar sgoldstino, 
the superpartner of the goldstino. The sgoldstino cou- 
plings are g g = M 3 /(2\/2F) and s 7 = M 77 /(2a/2F), 
where F is a parameter of supersymmetry breaking and 



Mr, 



Mi cos 9w+M 2 sin Oy/ with Mj the correspond- 



ing gaugino masses. Therefore, the sgoldstino coupling to 
photons is suppressed relative to gluons only by the "hi- 
erarchy among gauginos." Therefore, this is an example 
for a model where X couples to photons with a similar 
strength as to gluons. For M3 = 5M 77 = 500 GeV we 
obtain 



VF > 1.5 x 10 6 GeV 



instead of Eq. (Ef) and 



/ 10 20 eV 
V E x 



1/4 



M 



x 



Vf < 1.3 x 10 4 GeV 



Ex 



10 McV 



1/4 



10 20 eV 



(8) 



(9) 



instead of Eq. (0) . 

A variety of experimental limits on models with light 
sgoldstinos has been derived in 28 . In Fig. || we present 
the region of parameter space where sgoldstinos may act 
as UHECRs and are not excluded by other limits. This 
region corresponds to the upper region in Fig. [I]. 

If Ex = 10 21 eV or more, the allowed regions are 
larger, though no event of such energies has been ob- 
served. If g s = const/ A where A is the scale of new 
physics, then at const ~ 1 we have A = 10 2 -10 3 TeV. 
With E x = 10 11 GeV we have E cm = 300 TeV for in- 
teractions with protons. Certainly A should exceed this 
value if we want to use the nonrenormalizable interac- 
tions (|l|). For sgoldstinos we have M so ft ~ const F/A 
and A should be larger than E cm = 300 TeV. Note that 
F is a parameter of supersymmetry breaking and A is 
something like the scale of mediation of supersymmetry 
breaking which generally differs from y/F but should ex- 
ceed yF if const is of order 1. 



V. CONCLUSIONS 

We have suggested new (pseudo)scalar particles as Ul- 
tra High Energy Cosmic Rays beyond the GZK cutoff. 
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FIG. 2. Allowed region for the parameters (Mx, VF). The 
short-dashed line corresponds to the limit (Q), the long-dashed 
line to (Q) . Sgoldstinos with masses less than 10 keV (vertical 
solid line) are ruled out by the helium-burning life-time of 
horizontal-branch stars. 



Our analysis was particularly motivated by recent results 
suggesting that the sources of UHECRs are cosmologi- 
cally point-like Jjj and that at least some of the sources 
appear to be BL Lacertae jL4| at cosmological distances. 

We have calculated the required range of parameters 
characterizing these particles if we postulate that they 
should be produced in high-redshift sources, propagate 
through the Universe without decay or energy loss, and 
interact in the Earth's atmosphere strongly enough to 
produce extended air showers at energies beyond the 
GZK cutoff. The self-consistency of our analysis requires 
that the energy scale for new physics, which for SUSY 
models is the scale of mediation of supersymmetry break- 
ing, should be close to the UHECR center-of-mass energy 
with nucleons of E cnl = 300 TeV. 

As a specific example we studied light sgoldstinos. We 
considered restrictions on the parameters of the model 
which come from laboratory experiments and observa- 
tional data. We obtained the required region in parame- 
ter space of the model which obeys all existing limits. 

We note that our allowed region in Fig. 2 suggests 
that the supersymmetry breaking scale \/~F ~ 1-10 TeV. 
Hence our light sgoldstino model can be tested in searches 
for rare decays of JM and T and in reactor experiments 
(for details see Ref . |28j ) . This low scale of supersymme- 
try breaking may be also tested at new generation accel- 
erators like Tevatron and LHC. Also, sgoldstino contri- 
butions to FCNC and lepton flavor violation are strong 
enough to probe the supersymmetry breaking scale up to 
\fF ~ 10 4 TeV ^gj if off-diagonal entries in squark (slep- 
ton) mass matrices are close to the current limits in the 
MSSM. Thus our light-sgoldstino scenario for UHECRs 
allows only small flavor violation in the scalar sector of 
superpartners. 

Light (pseudo)scalars emerge not only in the context 
of supersymmetry, but also, for instance, in string the- 
ory and models with extra dimensions. Probably, such 
scalars also can serve as UHECRs if their effective cou- 
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pling with photons obeys the limits presented in sec- 
tion |I| 

Interpreting the UHECRs as new (pseudo)scalars is, 
of course, extremely speculative. However, we think it is 
noteworthy that such an interpretation is at all possible 
and self-consistent without violating existing limits. 
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